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AN  ADAPTIVE  MESH  TECHNIQUE  FOR  FINDING  THE 
BISTATIC  SPECULAR  POINT  ON  A  CONVEX  HULL 


INTRODUCTION 


One  of  the  most  computationally  expensive  routines  in  the  Naval  Undersea 
Warfare  Center’s  (NUWC’s)  High  Frequency  Target  Acoustic  Model  (HFTAM)  is  the 
routine  that  searches  for  a  bistatic  specular  point  along  a  hull-like  surface.  This  routine  is 
the  first  part  of  the  specular  scattering  calculation,  and  it  computationally  determines  the 
point  along  the  hull  surface  that  minimizes  the  distance  from  the  source  to  the  surface 
point  to  the  receiver.  This  point  will  be  unique  for  a  convex  shape  and  will  be  the  point 
whose  normal  bisects  the  vectors  from  the  point  to  the  source  and  from  the  point  to  the 
receiver.  This  bistatic  specular  point  is  used  by  other  routines  to  compute  the  specular 
scattering  from  the  hull-like  structure.  Monostatic  geometries  are  not  handled  separately 
and  are  only  viewed  as  a  case  where  the  source  and  receiver  happen  to  be  collocated. 

The  computational  routine  that  finds  the  specular  point  treats  the  hull-like  geometry 
description  with  the  source  and  receiver  locations  as  input  and  returns  the  coordinates  of 
the  specular  scattering  point  as  output.  For  the  remainder  of  this  report,  the  process  of 
finding  this  specular  point  for  a  fixed  source/receiver  pair  will  be  called  a  hull  search. 

This  report  begins  with  a  description  of  the  equations  used  to  model  convex  hulls, 
as  they  are  an  integral  component  of  the  solution  techniques.  Then,  a  detailed  explanation 
of  the  current  bisection  hull  search  technique  and  its  performance  is  presented.  This  is 
followed  by  a  description  of  the  new  adaptive  mesh  hull  search  technique  and  its  perfor¬ 
mance  relative  to  the  existing  technique.  Finally,  a  software  implementation  guide  is 
presented  for  incorporating  the  new  search  routine  into  the  HFTAM. 


GEOMETRY  OF  THE  CONVEX  HULLS 


The  surfaces  represented  as  convex  hulls  are  many  structures  that  can  be  internal 
or  external  to  a  target,  have  complex  curvature,  and  are  not  necessarily  represented  by 
canonical  shapes.  The  only  restriction  for  fitting  a  geometrical  shape  to  these  surfaces  is 
that  the  surface  must  be  representable  by  a  convex  shape  (or  a  portion  thereof).  Since  the 
majority  of  the  shapes  of  interest  are  submarine-like  (i.e.,  shaped  like  a  perturbed 
cylindrical  shell),  it  is  advantageous  to  represent  the  geometry  by  an  equation  that  is 
implicit  in  the  coordinates  and  is  of  the  form 


f(y,  ml(x),m2(x),...)  +  g(z,nx  (x),  n2  (x), . . . )  =  0. 


(1) 
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In  this  form,  the  functions  m;  and  tij  represent  parameters  that  vary  along  the  length  of  the 
hull  (hence  their  explicit  dependence  on  jc).  By  using  a  representation  of  this  form,  it  is  a 
direct  procedure  to  find  parameters  that  fit  a  given  geometry  to  the  equations  by  fitting 
“slices”  of  the  y-z  plane  at  fixed  values  of  jc  (i.e.,  identify  values  of  m/xi),  mfoi),  etc.,  for 
a  large  set  of  x  locations),  and  then  allowing  the  m;  and  parameters  to  vary  smoothly  as 
a  function  of  x.  In  this  context,  “smoothly”  implies  continuity  through  at  least  the  second 
derivatives  (C2  continuity). 

To  represent  the  types  of  curvature  often  seen  in  these  y-z  slices,  the  HFTAM  uses 
a  geometric  expression  for  each  slice  that  is  represented  by  a  hyperellipse:1 


y  \PW 
Kb(x) 


z-d(x)Yix) 


a(x) 


=  1, 


(2) 


where  the  parameter  d  allows  a  non-zero  offset  of  the  slice  in  the  z-direction.  An  example 
of  a  surface  defined  in  this  way  is  shown  in  figure  1.  The  expressions  for  the  parameters 
defining  this  particular  surface  can  be  found  in  appendix  A. 
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Figure  1.  Sample  Convex  Hull  Represented  by  the  Hyperellipttc  Formulation 


To  allow  for  a  greater  variety  of  surfaces,  the  hyperellipse  expression  can  be 
supplemented  by  a  set  of  cutouts.  These  are  regions  that  are  defined  in  x-z  space  and  are 
flagged  as  existing  on  the  mathematical  convex  hull  but  do  not  exist  on  the  real  surface. 
In  this  way,  partial  convex  hull  pieces  may  be  used  to  represent  complex  geometries. 
Hence,  the  specular  point  is  computed  as  if  it  were  a  complete  convex  hull,  and  then 
computational  logic  determines  whether  or  not  a  real  specular  point  exists  (which  occurs 
when  the  computed  point  is  not  in  the  flagged  region).  An  example  of  a  surface  defined 
with  a  cutout  is  shown  in  figure  2.  In  this  figure,  the  red  region  on  the  bottom  and  the 
blue  region  on  the  top  are  separated  by  a  cutout  shown  by  a  gray  slice. 

This  example  shows  how  a  piece  of  a  very  complex  curved  surface  can  be 
represented  by  a  convex  hull  with  associated  cutout  regions.  In  this  case,  the  red  section 
could  be  a  deck  structure  that  is  defined  separately  from  a  main  hull,  thus  allowing  a 
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nonconvex  geometry  to  be  represented  as  the  union  of  multiple  pieces  of  convex 
geometries.  Also,  the  hyperelliptic  exponent  p>2  for  the  red  section  (becoming  as  large 
as  20),  giving  the  convex  hull  section  a  squared-off  deck-like  appearance.  In  this  way, 
virtually  any  geometry  that  is  seen  in  an  undersea  structure  can  be  represented  by  a  series 
of  (potentially  partial)  convex  hull  sections.  The  expressions  for  the  parameters  that 
define  the  surface  and  the  cutout  region  shown  in  figure  2  are  found  in  appendix  B. 


y 

Figure  2.  Sample  Convex  Hull  with  a  Cutout  Region 


BISECTION  SEARCH  APPROACH 


The  current  technique  that  is  implemented  in  HFTAM  for  finding  the  specular 
point  was  provided  by  Lengua2,3  and  is  briefly  reviewed  here.  In  this  algorithm,  a  one¬ 
dimensional  bisection  search4  is  performed  in  the  x-direction.  At  each  step  of  the 
bisection  search,  the  source  and  receiver  points  are  projected  onto  the  y-z  plane.  An 
analytical  expression  for  the  minimization  of  the  distance  from  the  source  to  hull  to 
receiver  within  this  plane  is  then  evaluated. 

The  analytical  minimization  function  is  obtained  as  follows.  For  a  candidate  x 
location  given  by  x  =  xc ,  the  total  distance  from  the  source  to  a  point  on  the  hull  is  given 
by 

ps(y>  z)  =  Ws~  z?  +  (y,  -  yf  >  (3) 

where  the  subscript  s  corresponds  to  the  source’s  coordinates.  The  distance  to  the  receiver 
is  given  by  a  similar  expression  for  pr.  The  total  distance  is  minimized  for  a  given  x 
location  by  first  using  equation  (2)  for  the  fixed  x  value  to  solve  for  y  in  terms  of  z  as 
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With  this  definition  for  y,  the  minimization  problem  is  written  in  terms  of  the  single 
unknown  z  by  first  setting  the  distance  variation  to  zero  for  a  minimum  as 


4~(Ps+Pr)  =  0' 


(5) 


The  minimization  expression  in  equation  (5)  is  solved  by  combining  with  equation  (3)  to 
get 


Since  y  is  a  function  of  z  in  equation  (4),  this  expression  has  a  single  unknown  quantity  z. 
However,  the  expression  for  z  in  equation  (6)  will  be  transcendental  rather  than  a  simple 
polynomial  (due  to  both  y(z)  and  the  derivative  of  y  with  respect  to  z)  and  has  to  be 
solved  by  a  numerical  procedure.  The  HFTAM  program  currently  uses  a  bisection  search 
for  this  transcendental  evaluation. 

This  technique  is  mathematically  elegant  but  computationally  costly.  Further¬ 
more,  there  is  no  guarantee  of  convergence  of  the  algorithm,  because  the  convergence  is 
dependent  on  the  numerical  simplicity  of  the  analytical  minimization  problem.  If  the 
points  found  at  each  xc  vary  in  a  simple  manner  as  the  algorithm  moves  through  x  values, 
there  is  good  convergence.  However,  if  the  points  found  at  each  xc  vary  in  a  complicated 
way  (even  if  this  is  only  a  numerical  phenomenon),  the  convergence  may  be  poor  or  even 
nonexistent. 

The  bisection  search  technique  performs  well  when  the  hull  shape  is  nearly 
elliptical  (i .e.,p  =  2)  since  these  cases  tend  to  have  a  simple  analytical  expression  for  the 
minimum  distance  point  for  a  given  slice.  However,  for  the  cases  where  the  convex  hull 
geometry  is  used  to  represent  complex  shapes  other  than  hulls,  and  also  for  the  hulls  of 
many  diesel-electric  submarines,  the  hull  shapes  are  more  squared-off  at  the  comers 
( p  »  2)  rather  than  elliptical.  To  create  an  accurate  representation  of  these  nonelliptic 
geometries,  the  hyperellipse  formulation  uses  an  arbitrary  exponential  p(x),  which  may 
vary  along  the  length  of  the  hull.  It  is  in  these  cases  that  the  bisection  search  algorithm 
has  difficulty  in  converging  and  often  does  not  converge  in  the  maximum  number  of 
iterations  or  fails  to  converge  due  to  the  complexity  of  the  changes  in  the  parameter  p 
over  the  search  length.  Furthermore,  these  cases  tend  to  be  extremely  computationally 
intense  even  when  they  do  converge  to  a  solution.  This  computational  complexity  is 
caused  by  the  need  to  perform  a  secondary  bisection  search  at  each  candidate  x  value  as 
well  as  performing  a  large  number  of  noninteger  power  evaluations  in  the  analysis  of 
equation  (4)  and  its  derivative.  Thus,  there  are  costly  computations  that  occur  in  both  the 
evaluation  of  a  given  test  point  and  the  iteration  between  test  points.  It  should  be  noted. 


however,  that  the  computational  complexity  is  greatly  reduced  when  p  =  2.  In  that  case, 
the  expression  in  equation  (6)  can  be  rewritten  as  a  very  large  polynomial  equation  in  z, 
which  can  be  solved  quickly  and  accurately. 

To  illustrate  the  computational  complexity  of  these  routines  when  the  hull  shape 
is  nonelliptical  (i.e.,  p  *  2),  die  hull  search  evaluations  that  are  performed  for  two 
different  submarine-like  models  must  be  considered.  The  first  case  uses  a  model  called 
TestSub,  which  contains  both  hull  and  internal  hull-like  structures,  all  of  which  have  p 
identically  equal  to  two.  The  second  case  is  with  a  model  called  Omega,  which  uses  p 
values  that  not  only  are  different  from  two  but  also  vary  along  the  length  of  the  hull. 
Furthermore,  this  model  has  internal  hull-like  structures  that  have  p  values  other  than 
two.  To  perform  an  accurate  measure  of  the  computational  complexity  of  these  two 
models,  they  were  both  run  on  the  same  computer  (a  130-MHz  IBM  RS/6000)  over  the 
aspect  range  of -1 80°  to  +180°  at  a  0.00°  elevation  angle  in  a  monostatic  geometry.  The 
number  of  hull  search  calls  varies  as  a  function  of  aspect  as  various  internal  structures 
become  hidden  from  the  acoustic  path.  Thus,  the  total  time  spent  in  the  hull  search 
routine  over  all  aspect  runs  was  tallied  and  divided  by  the  total  number  of  hull  search 
calls  to  give  a  value  for  the  average  time  per  hull  search.  The  results  are  shown  in  table 
1.  It  is  clear  from  table  1  that  the  hull  search  routine  with  elliptical  curves  is  much  faster 
than  the  routine  for  nonelliptical  curves.  In  fact,  the  nonelliptical  routine  had  an  average 
search  time  that  was  over  12  times  longer  than  that  for  the  elliptical  hull.  Similar  results 
have  been  observed  for  other  models.  Thus,  the  bisection  hull  search  routine  has  been 
shown  to  be  very  slow  for  nonelliptical  hulls. 


Table  1.  Performance  of  the  Bisection  Hull  Search 


Model 

Number  of 
Hull  Searches 

Total  Time  in 
Hull  Searches 

Average  Time  per 
Hull  Search 

TestSub 

334 

0.20  sec 

0.60  msec 

Omega 

543 

4.15  sec 

7.64  msec 

For  real-time  (and  even  for  fast  wall-time)  simulation  purposes,  it  is  very 
important  to  make  all  of  the  routines  in  a  model  as  fast  as  possible  while  maintaining  a 
high  level  of  fidelity.  Furthermore,  a  guarantee  of  convergence  at  all  aspects  (where 
specular  points  exist)  is  a  requirement  for  simulations  because  the  model  is  used  only 
for  one  specific  aspect  at  a  time.  In  the  context  of  this  report,  this  is  referred  to  as  the 
robustness  of  the  numerical  algorithm.  The  robustness  of  the  bisection  search  algorithm 
was  addressed  by  examining  two  runs  of  the  Omega  model,  i.e.,  one  at  0.00°  elevation 
angle  and  a  second  at  0.01°.  The  results  of  these  two  runs  for  the  outer  hull  target 
strength  are  shown  in  figures  3  and  4.  It  is  obvious  from  these  plots  that  the  0.00° 
elevation  run  is  missing  the  beam  aspect  behavior  that  is  shown  in  the  0.01°  elevation 
run.  Some  analysis  of  the  runs  found  this  to  be  caused  by  a  lack  of  convergence  of  the 
hull  search  routine  in  the  former  run,  which  caused  the  target  strength  to  be  nonexistent. 
This  extreme  sensitivity  to  a  minor  change  in  elevation  angle  is  not  a  physical  effect  but 
merely  a  numerical  artifact  that  demonstrates  the  lack  of  robustness  of  the  bisection  hull 
search  technique. 
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Figure  3.  Target  Strength  for  Omega  Figure  4.  Target  Strength  for  Omega 

Model  Hull  with  Bisection  Hull  Search  Model  Hull  with  Bisection  Hull  Search 

at  0.00°  Elevation  Angle  at  0.01°  Elevation  Angle 


The  robustness  issue  is  a  primary  concern  only  in  simulation  applications;  it  is 
much  less  important  when  plots  of  behavior  over  a  range  of  aspects  is  the  required  output, 
such  as  for  target  strength  plots  or  “hourglass”  echo  plots.  However,  robustness  is  a 
critical  issue  when  the  output  of  the  model  is  going  to  be  used  as  input  for  a  processing 
scheme,  as  is  the  case  in  most  simulation  applications.  In  such  cases,  individual  runs 
(specific  aspect  angles)  are  performed  separately  and  the  individual  run  results  are  more 
important  than  the  composite.  A  numerical  search  technique  that  is  based  on  algebraic 
expressions  with  less  complexity  may  improve  this  robustness  issue.  When  such  an 
algorithm  is  developed  with  equivalent  or  better  speed  performance  than  the  bisection 
method,  it  will  be  a  practical  improvement. 


GEOMETRIC  MESH  SEARCH  APPROACH 


A  new  algorithm  for  the  bistatic  hull  search  problem  has  been  developed.  This 
algorithm  was  created  to  correct  the  robustness  and  speed  problems  of  the  bisection 
search  algorithm  for  simulation  purposes.  Robustness  issues  arise  when  the  geometry 
along  the  hull  shape  either  changes  extremely  rapidly  (sharp  comers)  or  very  slowly 
(cylindrical  sections)  and  the  bisection  algorithm  has  difficulty  converging.  Furthermore, 
the  speed  of  a  hull  search  routine  used  in  simulations  needs  to  be  relatively  insensitive  to 
the  geometric  complexity  of  the  specific  structure.  This  type  of  robustness  is  a 
requirement  of  the  user  community  to  ensure  that  the  computational  model’s  speed  is 
stable,  regardless  of  the  input  parameters  of  a  particular  run. 

To  handle  these  robustness  issues,  a  search  method  was  developed  based  on  local 
adaptive  refinement  of  geometric  meshes  along  the  convex  hull’s  surface.  The  advantage 
of  a  geometric  mesh  is  that  the  details  of  a  surface’s  complexity  lie  in  the  definition  of  the 
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mesh  points,  and  the  surface  between  the  mesh  points  is  represented  by  simple  curves. 
Thus,  once  the  mesh  is  generated,  all  of  the  complicated  equations  defining  the  surface 
have  been  replaced  by  a  table  of  mesh  points. 

For  the  acoustic  scattering  problem,  the  mesh  is  defined  to  a  level  of  detail  such 
that  the  local  radii  of  curvature  at  a  specular  point  are  determined  to  an  accuracy  of  at 
least  a  wavelength  squared  (X2).  This  requirement  is  driven  by  the  specular  scattering 
equation  from  the  Geometric  Theory  of  Diffraction5  and  its  numerical  evaluation.  This 
level  of  detail  can  be  guaranteed  only  if  the  space  between  nearest  mesh  points  is  no 
greater  than  the  wavelength  X.  Such  a  mesh  can  be  prohibitively  large;,  for  example, 
consider  a  cylinder  with  a  5-m  radius  and  a  length  of  50  m  being  ensonified  at  15  kHz. 
The  minimal  number  of  mesh  points  to  guarantee  a  one-wavelength  maximum  distance 
between  nearest  points  is  157,394  (this  does  not  even  account  for  the  endcaps).  Such  a 
large  number  of  mesh  points  is  not  practical  to  generate;  therefore,  a  meshing  scheme 
based  on  local  mesh  refinement  has  been  developed. 

Adaptive  mesh  refinement  is  practical  for  a  specular  point  searching  routine 
because  the  specular  points  are  a  local  phenomena,  and  the  mesh  only  needs  to  be  refined 
in  a  neighborhood  of  the  specular  point.  To  locally  determine  the  specular  point  at  a 
given  level  of  refinement,  the  mesh  is  evaluated  point- wise  and  the  minimum  combined 
distance  from  source  to  mesh  point  to  receiver  is  tabulated.  The  point  with  the  smallest 
combined  distance  is  the  closest  point  to  specular  in  the  given  level  of  refinement  for  the 
mesh  points  evaluated.  To  further  refine  the  mesh,  the  new  search  area  around  the  chosen 
point  must  be  made  large  enough  to  guarantee  convergence  as  the  mesh  refines.  In  a 
worst  case,  the  mesh  can  be  refined  by  at  least  a  level  of  one-fourth  (in  area)  at  each  step. 
In  this  scenario,  the  specular  point  will  be  known  to  within  an  accuracy  of 


after  n  steps  of  refining,  where  Ao  is  the  initial  convex  hull  surface  area,  SP  is  the  exact 
specular  point,  SP*  is  the  computed  specular  point,  and  dist( )  is  the  standard  I2  distance 
norm. 


The  locally  adaptive  meshing  technique  that  has  been  developed  for  the  hull 
search  problem  uses  a  mapping  of  the  convex  hull  surface  into  a  two-dimensional  space. 
The  mapping  was  specifically  derived  for  hyperelliptic  convex  hulls  represented  by 
equation  (2).  The  mapping  that  is  used  in  this  algorithm  translates  between  the  physical 
(x,y,z)  space  and  the  (x,  0)  space,  where  x  is  the  same  in  both  domains.  This  mapping  is 
given  by 


y(x,9)  =  ±b(x)- 


tantfp 

Iwj 


z{x,G)  =  d{x)±  a(x)  • 


1  \ 


where 
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\i!p(x) 


(9) 


P(0)  -  [l  +  |tan^|p(x) 
and  the  inverse  mapping  is  given  by 


6  =  arctan 


y  -  a{x) 


[z-d(x)\b(x)) 


(10) 


In  equation  (8),  the  sign  in  the  expression  fory  is  chosen  to  be  the  same  as  the  sign  of  6, 
and  the  sign  in  the  expression  for  z  is  chosen  to  be  negative  when  \6\  <  tt/2  and  positive 
otherwise.  In  all  of  these  expressions,  the  values  of  6*  are  assumed  to  be  in  the  range 
-n  <9<n . 


This  mapping  is  not  a  complete  geometric  mapping,  i.e.,  it  is  not  both  “one-to- 
one”  and  “onto”  mapping,  which,  in  general,  is  not  appropriate  for  numerical  meshing.6 
However,  the  mapping  is  onto  (i.e.,  every  point  in  one  domain  maps  to  the  other  domain 
and  vice  versa )  and  is  lacking  the  one-to-one  property.  This  point  may  seem  insignificant, 
but  it  is  important  since  there  are  multiple  points  in  the  (x,  6)  space  that  may  map  into  the 
same  point  in  the  (x,y,z)  space.  In  this  case,  it  is  not  a  problem  due  to  the  total  convexity 
of  the  shapes  being  modeled.  In  fact,  as  a  result  of  the  total  convexity,  it  is  guaranteed 
that  the  multiple  points  mapping  into  one  point  will  occur  at  exactly  two  locations  on  a 
convex  hull,  i.e.,  at  the  maximum  and  minimum  values  of  x,  where  the  hull’s  two- 
dimensional  representation  closes-up  on  itself.  Knowing  the  exact  location  of  these 
points  a  priori ,  a  correction  for  the  ambiguity  of  transforming  between  the  two-  and 
three-dimensional  spaces  at  these  points  is  made  with  simple  logic  in  the  computer 
program.  Furthermore,  since  the  mapping  is  onto,  there  is  no  problem  of  an  ill-defined 
point  in  either  space  (such  points  often  arise  due  to  numerical  inaccuracies)  that  would 
not  transform  to  the  other  space.  In  fact,  any  such  ill-defined  points  will  transform 
smoothly  between  the  two  spaces  (with  the  singular  exception  at  the  already  mentioned 
endpoints). 


In  the  adaptive  meshing  technique  that  has  been  developed,  an  initial  mesh  on  the 
convex  hull  is  performed  using  the  grid  spacing  of  20  points  in  a  rectangular  grid  in  the 
(x,  6)  space  that  is  given  by 


<i  = 


mm> 


(3*min  +  *max 


)  /  4,  (*rnin  -^max  )  ^  2, 


(*min  +  3xmax )  /  4, 


vmax 


0j={-7l/ 2,  0,  71 12,  7t  , 

where  xmm  is  the  minimum  value  of  the  defined  surface  in  the  x-direction  and  xmax  is  the 
maximum  value.  Each  of  these  mesh  points  is  transformed  to  the  full  three-dimensional 
space,  and  the  distance  from  the  source  to  the  mesh  point  dSy  and  the  receiver  to  the  mesh 
point  dRij  is  calculated  for  each  of  the  points  in  the  mesh.  The  point  whose  sum  of  these 
two  distances  represents  the  minimal  sum  is  returned  as  the  point  to  use  in  the  next 
iteration.  For  the  convex  shape  (the  only  shapes  under  consideration),  there  is  the 
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possibility  that  two  or  more  of  the  points  will  share  this  minimal  distance.  However,  this 
will  occur  only  for  neighboring  points,  which  is  another  benefit  of  the  convexity  of  the 
shapes  modeled.  For  this  reason,  the  iteration  scheme  includes  nearest-neighbor  points  in 
the  next  iteration. 

The  iteration  of  the  technique  is  performed  in  the  (x,  9)  space  with  transformations 
to  the  three-dimensional  space  performed  with  each  new  point  added  as  the  mesh  is 
refined.  These  transformations  to  the  physical  coordinate  system  are  the  most  computa¬ 
tionally  intense  portion  of  the  search  algorithm,  requiring  one  tangent  evaluation  and  two 
operations  of  raising  a  real  number  to  a  noninteger  power  for  each  mesh  point  (in  addition 
to  numerous  numerically  fast  operations).  Keeping  a  regularly  spaced  rectangular  mesh  in 
the  (x,  6)  space,  each  iteration  will  require  the  addition  of  16  new  mesh  points.  As  shown 
in  figure  5,  the  new  mesh  includes  the  entire  region  between  the  chosen  point  of  the 
previous  iteration  and  all  of  its  nearest  neighbors.  The  three-dimensional  representation  of 
each  of  the  new  points  is  used  to  compute  the  source  to  mesh  point  to  receiver  distances, 
which  are  compared  against  the  center  point.  Note  that  the  other  eight  points  in  the  refined 
mesh  patch  do  not  need  to  be  evaluated  since  they  were  used  at  the  previous  step  and 
determined  to  be  farther  than  the  center  point.  The  1 6  new  points  plus  the  center  point  that 
started  the  iteration  are  all  compared  for  minimum  distance  and  the  minimum  is  used  as  the 
center  point  for  the  next  iteration. 


Figure  5.  Diagram  of  the  Mesh  Adaptation  Procedure 


The  iterative  nature  of  the  hull  search  routine  allows  the  calculation  to  be  fairly 
quick,  but,  more  importantly,  its  speed  of  operation  can  be  predicted  without  assuming 
anything  about  the  convergence  of  any  particular  geometry.  In  practice,  for  a  convex  hull 
shape  with  length  L  and  maximum  diameter  D,  the  nominal  initial  area  can  be  assumed  to 
be  less  than 

Aq  =  kD{L  +  D),  (11) 
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which  is  the  area  of  a  cylinder  of  those  nominal  dimensions  with  spherical  endcaps. 
If  the  final  location  must  be  determined  to  within  a  wavelength  and  the  initial  area  is 
bounded  by  equation  (1 1),  then,  with  the  aid  of  equation  (7),  it  can  be  shown  that  the 
number  of  iterations  required  is 


log(rcD(Z  +  £>))-log^2A2  j 
log(4) 


(12) 


With  this  expression,  the  number  of  iterations  required  can  be  set  ahead  of  the  search  by 
knowing  only  the  maximum  diameter,  the  overall  convex  hull  length,  and  the  wavelength 
(or  frequency).  This  number  of  iterations  for  convergence  is  not  dependent  on  the 
complexity  of  the  shape  and,  therefore,  provides  an  improvement  in  terms  of  the 
computational  efficiency  of  the  hull  search  technique.  Also,  the  number  of  iterations 
required  increases  with  frequency,  but  only  as  log (f),  which  is  slower  than  linear  for  large 
frequencies. 

As  an  example  of  the  computational  requirements  of  adaptive  meshing,  consider 
again  the  case  of  a  cylinder  50  m  long  with  a  5-m  radius  being  ensonified  by  a  15-kHz 
pulse.  Direct  meshing  required  157,394  mesh  points,  not  including  the  endcaps.  For  the 
adaptive  mesh  technique,  equation  (12)  shows  that  nine  iterations  will  be  required.  Since 
20  mesh  points  are  generated  initially  and  an  additional  16  mesh  points  are  generated  at 
each  iteration,  this  problem  would  require  only  164  mesh  points  in  the  adaptive  technique. 
This  relatively  small  number  of  point  evaluations  (164  compared  to  more  than  150,000)  is 
what  makes  the  adaptive  meshing  technique  a  practical  approach.  Direct  meshing  would 
require  far  too  many  computations. 

To  assess  the  accuracy  and  robustness  of  the  adaptive  mesh  hull  search  method, 
comparison  runs  were  made  with  exactly  the  same  conditions  as  those  performed  for  the 
bisection  search  algorithm.  The  timing  comparison  is  shown  in  table  2  for  the  TestSub 
and  Omega  models. 

Table  2.  Performance  of  Adaptive  Mesh  Hull  Search 


Model 

Number  of 
Hull  Searches 

Total  Time  in  Hull 
Searches 

Average  Time  per 
Hull  Search 

TestSub 

334 

0.46  sec 

1.38  msec 

Omega 

543 

0.85  sec 

1.56  msec 

In  contrast  to  the  bisection  hull  search,  the  adaptive  mesh  hull  search  technique 
has  a  much  more  consistent  average  time  per  hull  search.  With  this  technique,  the 
Omega  model  average  time  per  hull  search  is  only  13  percent  longer  than  TestSub’s; 
whereas  for  the  bisection  search,  the  difference  was  1173  percent.  However,  the  average 
time  per  hull  search  for  the  elliptical  case  of  TestSub  (where  the  bisection  technique 
converges  rapidly)  is  more  than  twice  as  long  with  the  adaptive  mesh  hull  search. 
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Two  more  models  were  run  with  both  hull  search  techniques  to  confirm  that  the 
timing  conclusions  were  consistent  and  not  based  on  any  peculiarities  of  a  specific  model. 
The  two  models  added  were  called  Yellow  (a  model  with  all  elliptical  hull  shapes)  and 
XYZ  (a  model  with  many  hyperelliptical  hull  shapes).  The  results  of  the  timing  runs  for 
all  four  models  with  both  hull  search  techniques  are  summarized  in  table  3,  which 
presents  only  the  average  time  per  hull  search. 


Table  3.  Performance  of  Bisection  and  Adaptive  Mesh  Hull  Searches 


Model 

Elliptical 

Number  of 
Searches 

Bisection 
Average  Time 

Adaptive  Mesh 
Average  Time 

TestSub 

Yes 

334 

0.60  sec 

1.38  msec 

Yellow 

Yes 

297 

0.68  sec 

1.31  msec 

Omega 

No 

543 

7.64  sec 

1.56  msec 

XYZ 

No 

629 

4.75  sec 

1.86  msec 

From  the  results  shown  in  table  3,  it  is  evident  that  the  adaptive  mesh  search 
algorithm  is  much  more  consistent  in  runtime  than  the  bisection  search.  Furthermore,  the 
worst  case  of  the  adaptive  mesh  hull  search  (1 .86  msec)  is  much  faster  than  the  worst 
case  of  the  bisection  hull  search  (7.64  msec).  This  is  because  the  bisection  hull  search  is 
much  slower  for  nonelliptical  hulls,  whereas  its  performance  for  elliptical  hulls  is  very 
fast.  By  contrast,  the  adaptive  mesh  hull  search  is  only  slightly  slower  for  nonelliptical 
hulls  than  it  is  for  elliptical  hulls.  This  slowdown  is  attributed  to  the  increased  complexity 
of  point  transformations  from  two-dimensional  to  three-dimensional  space  when  the 
exponents  are  not  equal  to  two  (see  equation  (8)). 

The  consistency  in  average  hull  search  time  that  is  found  in  the  adaptive  mesh 
search  algorithm  is  one  of  its  improvements  over  the  bisection  search.  It  should  be  noted 
that  the  increase  in  computational  time  using  the  adaptive  mesh  for  elliptical  hulls  (over 
the  bisection  search)  is  not  large  enough  to  outweigh  the  performance  gains  found  with 
nonelliptical  hull  sections. 

The  second  performance  improvement  that  is  gained  by  using  the  adaptive  mesh 
search  algorithm  is  in  the  robustness  of  the  computation.  Running  the  same  two  cases 
that  were  shown  in  figures  3  and  4  (Omega  model  outer  hull  with  bisection  at  0.00°  and 
0.01°  elevation  angles,  respectively)  showed  that  the  adaptive  hull  search  algorithm 
produced  a  consistent  result  over  the  small  change  in  elevation  angle.  The  results  of  hull- 
only  target  strength  runs  for  these  cases  are  shown  in  figures  6  and  7.  From  the  plots,  it  is 
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Figure  6.  Target  Strength  for  Omega  Figure  7.  Target  Strength  for  Omega 

Model  Hull  with  Adaptive  Mesh  Hull  Model  Hull  with  Adaptive  Mesh  Hull 

Search  at  0. 00°  Elevation  Angle  Search  at  0. 01  °  Elevation  Angle 


obvious  that  the  adaptive  mesh  search  routine  provides  a  robust  result  that  does  not  vary 
drastically  with  extremely  small  angle  changes. 

The  peculiar  behavior  of  the  target  strength  plot  near  the  0°  aspect  angle  is  a 
function  of  the  squared-off  nature  of  the  bow.  This  type  of  behavior  can  be  seen  only  in 
complex  hull  shapes  and  is  not  seen  in  simple  canonical  hulls.  This  feature  is  also  not 
captured  in  the  bisection  search  algorithm  since  the  bow  is  at  the  endpoint  of  the  two- 
dimensional  coordinate  surface  and  the  bisection  search  algorithm  requires  many  more 
iterations  than  the  maximum  allowed  to  converge  there  for  nontrivial  hull  shapes.  The 
behavior  of  the  target  strength  plot  in  the  stem  section  (beyond  110°  aspect  angle)  is 
something  that  is  peculiar  to  a  complete  convex  hull.  This  effect  is  from  the  small  cap¬ 
like  section  at  the  stem  required  to  make  the  shape  complete.  If  the  effect  is  not  desired,  it 
can  be  removed  by  adding  a  cutout  to  the  hull  geometry  definition  in  that  region. 

The  accuracy  of  an  overall  acoustic  scattering  model  such  as  the  HFTAM  will  be 
affected  by  a  change  in  the  hull  searching  algorithm.  Since  the  HFTAM  has  already  been 
validated  at  the  peak  target  strength  level,  the  impact  of  any  changes  at  the  peak  target 
strength  level  needs  to  be  documented.  Peak  target  strength  calculations  for  the  entire 
structure  (outer  hull,  sail,  internals,  etc.)  were  computed  for  the  HFTAM  with  the 
bisection  hull  search  routine  and  with  the  adaptive  mesh  hull  search  routine  for  all  four 
targets  studied.  The  results  of  the  peak  target  strength  comparisons  can  be  seen  in  figures 
8  through  15.  All  of  these  results  were  plotted  in  a  polar  format  with  a  radial  direction 
displaying  level  in  decibels  and  the  rotational  direction  representing  aspect  angle  (off- 
bow)  in  degrees. 
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As  can  be  seen  in  figures  8  through  1 1,  the  two  models  that  had  all  elliptical  hull 
cross  sections  (TestSub  and  Yellow)  produced  similar  peak  target  strength  results  for  the 
two  different  hull  search  techniques.  This  is  expected  because  the  bisection  algorithm 
has  no  problems  with  purely  elliptical  hull  geometries.  The  adaptive  mesh  hull  search 
technique  always  converges  to  a  point  that  is  near  the  true  specular  point,  so  it  is  believed 
to  be  a  more  accurate  technique.  However,  the  existing  HFTAM  with  the  bisection 
search  algorithm  has  also  been  validated  against  acoustic  data  and,  therefore,  must  be 


Figure  8.  Peak  Target  Strength  of  TestSub 
Model  Using  Bisection  Hull 
Search  Technique 


Figure  9.  Peak  Target  Strength  of  TestSub 
Model  Using  Adaptive  Mesh  Hull 
Search  Technique 
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Figure  10.  Peak  Target  Strength  of  Yellow  Figure  11.  Peak  Target  Strength  of  Yellow 
Model  Using  Bisection  Hull  Model  Using  Adaptive  Mesh  Hull 

Search  Technique  Search  Technique 
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assumed  to  be  accurate  for  the  shapes  that  have  been  used  for  the  model  validation.  The 
existing  validations  were  performed  with  target  shapes  that  are  purely  elliptical,  so  a 
comparison  of  figures  8  and  9,  as  well  as  figures  1 0  and  1 1 ,  shows  that  existing  validation 
data  are  still  applicable  to  the  adaptive  mesh  hull  search  algorithm.  Thus,  the  adaptive 
mesh  hull  search  algorithm  is  at  least  as  valid  as  the  bisection  algorithm,  when  compared 
with  previous  validations. 

For  nonelliptical  hull  geometries,  such  as  those  found  in  the  Omega  and  XYZ 
models,  the  adaptive  mesh  algorithm  computes  a  specular  point  for  all  cases.  However, 
there  are  times  when  the  bisection  algorithm  fails  to  converge  in  the  Omega  and  XYZ 
models  and,  therefore,  returns  no  target  strength.  Thus,  the  total  net  peak  target  strength 
for  all  highlights  may,  for  certain  aspects,  be  lower  in  the  bisection  algorithm  than  in  the 
adaptive  mesh  algorithm.  The  plots  in  figures  12  through  15  confirm  this.  In  figures  12 
and  13,  for  the  Omega  model,  the  difference  in  peak  target  strength  is  very  small  between 
the  two  techniques.  This  is  due  to  the  hull  not  being  the  dominant  highlight  over  most 
aspect  angles  (the  exception  being  near  bow  and  beam).  In  figures  14  and  15,  for  the 
XYZ  model,  there  is  a  considerable  difference  in  peak  target  strength  between  the  two 
techniques,  which  can  be  attributed  to  the  dominance  of  hull  highlights  for  this  particular 
target.  Also,  the  behavior  near  the  stem  (180°)  that  is  exhibited  in  figure  1 5  is  caused  by 
computing  the  reflection  near  the  tail  section  of  the  outer  hull,  a  region  where  the 
bisection  algorithm  fails  to  converge  (see  figure  14). 

The  discrepancies  found  between  the  two  techniques  in  figures  14  and  15  should 
not  affect  any  existing  validations  of  the  HFTAM  since  these  are  based  on  the  XYZ 
model.  The  XYZ  model  has  not  yet  been  validated  at  the  peak  target  strength  level 
against  acoustic  data;  therefore,  it  is  not  imperative  to  make  the  two  techniques  agree  at 


Figure  12.  Peak  Target  Strength  of  Omega  Figure  13.  Peak  Target  Strength  of  Omega 
Model  Using  Bisection  Hull  Model  Using  Adaptive  Mesh  Hull 
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Figure  14.  Peak  Target  Strength  ofXYZ 
Model  Using  Bisection  Hull 
Search  Technique 
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Figure  15.  Peak  Target  Strength  of  XYZ 
Model  Using  Adaptive  Mesh 
Hull  Search  Technique 


the  peak  target  strength  level.  It  is  anticipated  that  the  results  with  the  adaptive  mesh  hull 
search  technique  (figure  15)  will  be  more  accurate  than  the  results  for  the  bisection  hull 
search  technique  (figure  14). 

The  peak  target  strength  validation  metric  is  not  the  only  validation  metric  avail¬ 
able,  but  it  has  been  presented  here  to  show  the  validity  of  the  new  hull  search  technique. 
This  analysis  provides  a  connection  with  the  existing  validation  databases.  The  new 
technique  should  also  be  validated  by  specific  users  with  regard  to  their  needs,  such  as 
highlight  location,  time  echo  characteristics,  etc.  The  comparisons  shown  in  this  report 
are  illustrative  only  and  do  not  represent  a  complete  validation. 


SOFTWARE  CHANGES  FOR  IMPLEMENTATION 


To  implement  the  adaptive  mesh  specular  point  search  algorithm,  a  number  of 
changes  needed  to  be  made  to  the  HFTAM  software.  These  changes  have  all  been  made 
and  tested  on  HFTAM  (V ersion  1 .5)  to  produce  the  results  shown  in  this  report.  It  is 
expected  that  the  adaptive  mesh  hull  search  routine  will  become  the  default  hull  search 
technique  in  the  next  version  of  the  HFTAM  (V  ersion  1.6). 

The  HFTAM  software  changes  were  all  written  in  ANSI  C  and  conform  to  the 
programming  style  used  in  the  other  HFTAM  routines.  (All  of  the  software  changes  that 
are  documented  in  this  report  are  maintained  by  the  author  of  this  report.)  The  programs 
were  tested  by  creating  a  new  model  routine  bisim2.c  that  was  very  similar  to  the  existing 
bisim.c  with  the  exception  of  calling  up  the  new  hull  search  in  place  of  the  old  hull 
search.  Table  4  lists  the  programs  that  were  created  for  the  new  hull  search  routine. 
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Table  4.  Additional  HFTAM  Software  Routines  Written 
for  the  Adaptive  Mesh  Implementation 


Name 

Description 

HullSeach2.c 

Main  program  for  adaptive  meshing  used  in  the 
same  way  as  HullSearch.c. 

PatchEval.c 

Finds  the  point  with  minimal  distance  for  a 
given  level  of  mesh  refinement. 

Xtheta2XYZ.c 

Computes  the  three-dimensional  coordinates  of 
a  point  from  its  two-dimensional  representation. 

ddelli2.c 

Computes  parameter  values  only  (without 
derivatives)  for  a  fairing  ellipse  section. 

ddpoly2.c 

Computes  parameter  values  only  (without 
derivatives)  for  a  polynomial  section. 

GetFitParameters2  .c 

Finds  ( a,b,d,p )  for  a  given  x  value. 

Table  5  lists  the  routines  that  required  modification  to  accommodate  the  new  hull 
search  routine. 


Table  5.  HFTAM  Software  Routines  That  Are  Changed 
for  the  Adaptive  Mesh  Implementation 


Name 

Description 

Echoes.c,  bisim.c 

Recompiled  to  be  linked  with  new  routines. 

m_protos.h 

Modified  to  include  the  input/output  and  call- 
list  structure  of  each  of  the  new  routines  shown 
in  table  4. 

Hull.c 

Changed  to  call  HullSearch2  instead  of  calling 
HullSearch. 

When  the  new  hull  search  routines  are  in  place,  some  old  search  routines  will 
become  obsolete.  Routines  that  were  specific  to  the  bisection  hull  search  routine  and  that 
will  no  longer  be  required  in  an  implementation  of  HFTAM  are  listed  in  table  6. 
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Table  6.  Software  Routines  That  Become  Superseded  in  the 
Adaptive  Mesh  Implementation 


Name 

Description 

GetCandidatePoint.c 

Finds  iyf)  coordinates  for  a  sample  x  value. 

HullSearchx 

Replaced  by  HullSearch2.c. 

YZEqn.c 

Performs  the  bisection  search  on  a  fixed  x  slice 
for  the  analytical  distance  minimization 
problem. 

CONCLUSIONS 


A  new  adaptive  mesh  technique  for  searching  for  a  bistatic  specular  point  on  a 
convex  hull  surface  has  been  developed  for  the  NUWC  HFTAM.  The  existing  method  of 
bistatic  specular  point  searching  in  HFTAM  uses  a  bisection  search  with  numerical 
evaluations  of  a  transcendental  equation  at  each  bisection  step.  This  method  has  been 
shown  to  lack  computational  robustness. 

The  adaptive  mesh  technique  has  been  implemented  with  the  HFTAM  and  has 
been  shown  to  be  more  robust  than  the  bisection  technique.  Furthermore,  the  speed  of 
computation  for  the  new  technique  is  independent  of  the  complexity  of  the  geometric 
shape  being  modeled.  The  bisection  technique,  in  contrast,  was  shown  to  perform  very 
slowly  when  nonelliptical  hull  sections  were  used.  For  many  of  the  complex  shapes  that 
are  defined  to  a  high  level  of  geometric  fidelity,  the  nonelliptical  representations  are 
necessary  and,  thus,  there  is  a  need  for  the  new  hull  search  technique. 

The  adaptive  mesh  hull  search  technique  is  expected  to  become  the  hull  search 
used  in  a  future  release  of  the  HFTAM.  It  has  been  shown  to  improve  robustness  of  the 
current  technique  and  to  improve  the  speed  of  the  search  process  for  complex  geometries. 
The  impact  of  the  new  search  routine  on  the  HFTAM  has  been  shown  in  terms  of  overall 
peak  target  strength  as  well  as  some  individual  hull  target  strength  comparisons.  The 
model  shows  results  that  are  at  least  as  accurate  as  the  existing  technique.  These 
improvements  in  the  hull  search  routine  will  improve  the  speed  and  fidelity  of  the 
HFTAM  for  future  simulation  studies. 
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APPENDIX  A 

HYPERELLIPSE  PARAMETERS  FOR 
SAMPLE  CONVEX  HULL  OF  FIGURE  1 

This  appendix  contains  the  equations  defining  the  hyperellipse  parameters  for  the 
surface  shown  in  figure  1.  The  convex  surface  of  this  shape  is  defined  on  the  region 
-  55.1688  <  x  <  55.1688. .  The  equations  defining  the  surface  are  the  generalized 
hyperellipse  equations  in  the  form  of  equation  (2),  and  the  parameter  a  is  given  by  the 
following: 

18.7558  +  0.47683*  +  0.021646*2 


+(4.2292  x  10"4)*3  +  (2.4481  x  10“6)*4 

for 

-55.1688  <*<-31.5, 

-  3 .9304 1  +  V-  0.02805*2  -  0.66686*  +  76.3 1 1 07 

for 

-31.5  <*<-11.8872, 

5.0292 

for 

-11.8872  <*<36.8808, 

2.9442  +  V-  0.027338*2  +  2.01651*  -  32.8381 

for 

36.8808  <  *  <  44.5069, 

-  22.027765  +  1.32234*  -  0.016266*2 

for 

44.5069  <  *  <  49.9933, 

V-  0.35705*2  + 1 .48663*  +  26.65693 

for 

49.9933  <*<55.1688. 

For  this  particular  shape,  b  =  a  for  all  * .  Furthermore,  for  this  shape,  the  parameters  d 
and  p  satisfy  the  following: 

d  =  0,  for  -55.1688  <*  <55.1688 
p  =  2,  for -55.1688  <*<55.1688. 

Therefore,  all  of  the  cross  sections  of  the  shape  are  circles  (since  a  =  b  and  p  =  2)  that 
are  centered  about  the  *-axis  (since  d  =  0 ). 


A-l/(A-2  blank) 


APPENDIX  B 

HYPERELLIPSE  PARAMETERS  FOR 
SAMPLE  CONVEX  HULL  OF  FIGURE  2 


This  appendix  contains  the  equations  defining  the  hyperellipse  parameters  for  the 
surface  shown  in  figure  2.  The  convex  surface  of  this  shape  is  defined  on  the  region 
-  9.89  <  x  <  1 1 .0.  The  equations  defining  the  surface  are  the  generalized  hyperellipse 
equations  in  the  form  of  equation  (2).  The  parameter  d  for  this  shape  is  varying,  and 
therefore,  the  centerline  of  the  cross  sections  is  nonconstant  throughout  the  x-range.  In 
this  case,  the  d  parameter  takes  on  the  value 


-3.01835  +  (3.7245 x  10‘3)x  - (7.5864 x  10_4)x2  for  -9.89  <x< -2.2307, 
-2.8479  +  0.53 107*  +  0.62639x2 

+  0.3 5628x3+0.09839x4  +  0.01 064x5  for  -2.2307  <x  <-1.3135, 


-3.026132 


for  -1.3135  <  x  <11.0. 


Using  this  variable  d  parameter,  the  cross  sections  are  defined  in  two  components.  The 
first  components  are  labeled  with  a  subscript  1  (i.e.,  a\,  b\  and  p\)  and  are  valid  for  points 
where  z  >  d(x) .  The  second  components  are  labeled  with  a  subscript  2  (i.e.,  b%,  and 

P2)  and  are  valid  for  points  where  z  <  d(x) .  To  have  a  smooth  (up  to  C2  smoothness  is 
guaranteed)  junction  between  these  two  parts  of  each  cross  section,  the  b  parameters  are 
forced  to  be  the  same,  bx  =  b2. 

For  the  convex  hull  shape  in  figure  2,  the  cutout  region  will  be  defined  to  contain 
the  entire  region  of  the  subscript-1  parameters  (where  z  >  d(x)),  as  well  as  a  portion  of 
the  subscript-2  parameters.  Thus,  the  values  of  the  subscript- 1  parameters  are  somewhat 
arbitrary.  For  convenience,  the  values  of  the  a  parameters  can  be  forced  to  be  the  same 
by  making  a\  the  same  as  ai.  Also,  for  further  simplicity,  the  value  of  p\  is  assumed  to 
be  constant  at  p\ -2  throughout  the  entire  x-domain.  The  parameter  a\  (which  for  this 
shape  is  the  same  as  ai)  is  given  by  the  following  expression: 


B-l 


-  72.5 1 13  - 1 5.8452*  -  0.86036*2 

for 

-9.89  <*<-9.75, 

-  333247.6  - 176750.3*  -  37483.41*2  -  3972.92*3 
-210.459*4  -4.4575* 5 

for 

-9.75  <*<-9.175, 

1 .38480  -  0.054875*  -  0.01 71 89*2 

for 

-9.175  <*<-2.23074, 

1 .96 1 1 9  + 1 .6046*  + 1 .9074*2  +1.1 0659*3 
‘  +  0.3 1345*4 +0.034934*5 

for 

-  2.23074  <*  <-1.3135, 

1.43304 

for 

-1.3135  <*<9.4, 

50240.1  -  24949.4*  +  4948. 6*2  -  490.03 8*3 
+  24.2284*4 -0.47851*5 

for 

9.4  <  *  <  10.7, 

-  6.8818  +  2.3587*  -  0.15923*2 

for 

10.7  <  *  <  11.0. 

The  last  parameter  that  remains  to  be  defined  for  this  convex  hull  shape  is  p\.  The 
parameter  p\  is  defined  by  the  following  expression: 


22. 1 338  +  0.94835*  -  0.3 1 932*2 
-0.040468*3  -(l  .93717  x  10-3  )c4 

for 

-9.89  <*<-2.2307, 

23.3906  + 1 1.8417*  + 19.1948*2  + 14.2006*3 
+  4.67320*4+  0.57471*5 

for 

-2.2307  <*<-0.8135, 

20.657 

for 

-0.8135  <*<5.3736, 

40279.12  -  33263.6*  + 10932.12*2  - 1786.003*3 
+  145.033*4  -4.68386*5 

for 

5.3736  <*<6.3718, 

784.1 19  -  340.423*  +  56.1141*2 
-4.11922*3 +0.1 13 12*4 

for 

6.3718  <  *  <  1 1.0. 

The  preceding  expressions  define  a  complete  convex  hull  shape.  To  obtain  the 
shape  found  in  red  in  figure  2,  the  cutout  region  needs  to  be  defined.  The  cutout  region 
for  this  convex  hull  shape  is  defined  by  all  values  of  z  that  are  greater  than  the  following: 


B-2 


-3.02819 

for 

-5.36825  +  0.10819* 

for 

0.001  l  +  00154x  +  0.1448x2 

-0.0115X3  +  0.00022x4 

for 

-6.7875  + 0.02 169x 

for 

-0.2265 

for 

9.73777  <x<  21.63, 
21.63  <x  <27.35, 

27.35  <x<  29.75, 
29.75  <x  <30.25, 

30.25  <x£  30.95. 


B-3/(B-4  blank) 
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